Tuning MgB 2 (0001) surface states 
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Surface state localization and hybridization on B, Mg, and Li terminated MgB 2 (0001) surfaces are studied 
within density functional theory (DFT). For the B terminated surface we find a low energy B o~\ surface state, 
a sp z surface state, and B tr 2 and <r 3 quantum well states, which are 90% localized in the topmost B layer. 
Our results demonstrate that by charging the B atomic layer, either by changing the surface termination or 
through electro-chemical doping, the B cr surface states are shifted down in energy, filled, and delocalize through 
hybridization as they cross the bulk MgB 2 bulk bands. On the other hand, the sp- surface state should be shifted 
up in energy, emptied, and gain an increasingly metallic s character by adding a Mg, Mg +1 , or Li terminating 
atomic layer. These results clearly show both the robust nature of MgB 2 (0001) surface states, and how their 
localization and energy range may be tuned by surface termination and charging, with implications for the 
superconducting and plasmonic behaviour of MgB 2 . 

PACS numbers: 73.20.-r, 74.25 Jb, 74.70.Ad 



I. INTRODUCTION 

Since 2001 when superconductivity in MgB2 was first 
discovered 1 , many calculations of its bulk and surface elec- 
tronic structure have been performed. In general, these calcu- 
lations showed a good agreement with ARPES experiments, 
suggesting that MgB2 is not a strongly correlated system, 
and that superconductivity is probably a consequence of some 
strong electron-phonon coupling mechanism. In other words, 
MgB2 appears to be an extreme case of an Eliashberg super- 
conductor. Indeed, many studies show strong coupling of cr 
electrons and e^ g phonons in the B layer, and a wide super- 
conducting gap in the cr band has been measured 2 . Also, since 
the recent fabrication of pristine ultra-thin MgB2 films 3 , sur- 
face enhanced superconductivity has become an increasingly 
active area of research 3 . 

These theoretical investigations greatly contributed to the 
understanding of the peculiar bulk/surface electronic struc- 
ture of the seemingly simple MgB2 crystal. Bulk electronic 
structure calculations 4-7 show that MgB2 has five main bands, 
namely completely filled cr\ B, partially filled <t 2 , <t 3 and p z 
B, along with completely unoccupied s Mg bands. Bonding 
within the B layers is mostly covalent {sp 2 hybridization) and 
ionic/metallic between B-Mg-B layers. 

Due to the Mg between the B planes, the nature of in-plane 
covalent bonding is not typical, as in graphite for example. 
Specifically, the presence of Mg enhances the inter-layer over- 
lap and induces strong hybridization between dangling B p z 
and Mg s orbitals (sp z bond). This causes an upward shift of 
the Mg s and downward shift of the B p, band. In other words, 
this induces an emptying of the Mg s and filling of the B p z 
bands. In this way Mg atoms donate electrons interstitially 
(between B-Mg layers) and become positively ionized. This 
causes an additional downward shift of the B n band relative 
to the B cr bands and charge transfer from the cr to the n band. 

In this way MgB2 has a small hole doping of the higher en- 



ergy cr 2 and o"3 bands (at the F point) and bonds within the B 
layers become mixed covalent/metallic. Such hole doping in 
the B layers enhances the superconductivity. This is similar to 
the case of intercalated graphites 8 , where doping (inter-plane 
intercalation) by alkali metals causes the graphite to become 
superconducting. 

The above description is also consistent with charge den- 
sity distributions calculated in Ref. 4, where it was seen that 
Mg atoms are strongly ionized, but electrons are donated in- 
terstitially, rather than directly to the B layers, and inter-layer 
bonding is much more metallic than ionic. Also, it has been 
shown that less charge is participating in the in-plane B-B 
bonding in MgB2 compared to the C-C bonding in primitive 
graphite. This suggests that Mg induces a covalent bonding 
in the B layer which is not completely saturated and is weak- 
ened, i.e. slightly metallic. 

On the other hand, the electronic structure of the MgB2 sur- 
face is much more delicate. Due to the weak dispersion of 
bulk B cr bands in the (0001) direction, there exist two wide 
gaps in the projected bulk band structure of MgB2 at the F 
point. This allows for the formation of many different types 
of surface states. Indeed, several studies of the MgB2(0001) 
surface 9-12 showed the existence of various surface and sub- 
surface states, whose positions and existence depend on crys- 
tal termination and coverage 13 . 

For the case of a B terminated surface, every bulk band has 
its own surface state band which follows the upper edge of the 
corresponding projected bulk band. All these states are mostly 
localized in the first few B layers. For example, the surface 
state band for B terminated MgB2 reported in Ref. 1 1 has p z 
symmetry and is located on the upper edge of the projected p z 
bulk band. It's energy at the T point in the surface Brillouin 
zone (SPZ) is about -2.7 eV and its charge density is mostly 
distributed throughout the first three B layers. 

In the case of the Mg terminated surface, there is one sur- 
face state whose energy at the F point is about -2 eV. This 
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surface state has sp z symmetry and is mostly localized above 
the surface plane in the vacuum region (supra-surface). 

However, ARPES measurements which followed 14 
showed a good agreement with bulk electronic structure 
calculations 5 7 , although the peak at 0.5 eV in the normal 
photo-emission does not agree with any of the above men- 
tioned surface states. The calculations of Servedio et al. 12 
showed good agreement but are not particularly relevant, as 
the surface potential barrier was there modeled by a step 
potential, so that the binding energy of the surface state was 
strongly dependent on the position of the potential step. On 
the other hand, calculations by Profeta et al. 13 showed that 
the sp z surface state at the F point strongly depends on the 
termination and coverage of the MgB2 crystal, indicating that 
in the experiment 14 the surface termination may have been 
mixed or had significant surface contamination. Uchiyama 
et al. 2 paid much more attention to surface preparation, 
and in their measurement the normal photo-emission peak 
appeared at -1.0 eV. In Refs. 15-17 an 18 monolayer (ML) 
thick (18 ML B and 18 ML Mg layers) MgB 2 film was 
deposited an a Mg(0001) substrate using the Molecular Beam 
Epitaxy (MBE) co-deposition technique. Following this it 
was determined (using several techniques including XPS, 
LEED+XRD, X ray absorption) that indeed, there was no 
surface contamination, and ARPES measurements were then 
performed. The results showed an excellent agreement of 
the projected bulk band structure with bulk calculations 5 ' 6 . 
Agreement with surface band structure calculations 11 was 
also quite good. For example, peaks at about -1.6 eV and 
-3.2 eV are measured which correspond to sp z surface states 
in Mg and B terminated surfaces respectively. Measurements 
in Refs. 5,6 also suggested that the most stable surface 
termination is Mg, as was anticipated in Refs. 10,13. 

Recent experimental realizations of the MgB 2 surface have 
been brought to an increasingly high level of precision 2 ' 3 ' 15-17 . 
These experiments were aimed at investigating the fine super- 
conductivity properties of MgB 2 , e.g., multiple superconduc- 
tivity gaps 2 , and surface enhanced superconductivity 3 . How- 
ever, such clean MgB 2 surfaces, with a very well defined elec- 
tronic structure, can also be an excellent starting point for 
the investigation of both optical and dielectric MgB 2 surface 
properties. It also should be noted that in Ref. 3 ultrathin 
MgB 2 films down to 7.5 nm were fabricated, which is quite 
close to the thicknesses used in the calculations shown herein. 

Interplay between 2D and 3D electronic structures give 
rise to the unique dielectric 18-21 and extraordinary optical 
properties 22 ' 23 of the MgB 2 crystal. On the other hand, be- 
cause of the plethora of surface, localized and subsurface 
states, the MgB 2 surface dielectric and optical properties may 
be even more interesting. Before providing such theoretical 
calculations we require an accurate description of the surface 
electronic structure. 

For this reason, we will focus here on the investigation 
of three (B, Mg and Li) terminated MgB 2 surface electronic 
structures by comparing the results of three different (two 
plane wave and real space) ab initio density functional the- 
ory (DFT) methods. In so doing, we may demonstrate how 
by changing the surface termination, charging, and doping the 



MgB 2 (0001) we may tune both the localization and energy 
range of the surface and subsurface states, with potential ap- 
plications in the areas of plasmonics and superconductivity. 

The paper is organized as follows. In Sec. II we describe the 
ab initio methods used for the surface electronic structure cal- 
culations, the description of surface formation energies, and 
develop criteria to distinguish between surface, subsurface, 
vacuum, and bulk states. The numerical results are presented 
in Sec. Ill for B, Mg, and Li terminated surfaces, with results 
for each termination in separate subsections. These results are 
then discussed in more detail in Sec. V, where we show how 
surface termination and electro-chemical doping may be used 
to tune MgB 2 (0001) surface states. This is followed by a con- 
cluding section. 



II. METHODOLOGY 

Structural optimization and electronic structure calcula- 
tions for each MgB 2 (0001) surface have been performed by 
combining three different ab initio codes. Specifically, we 
compare our own pseudopotential based plane wave DFT 
code (SPPW) and the plane wave self-consistent field DFT 
code (PWscf) belonging to the Quantum Espresso (QE) 
package 24 , with a real-space projector augmented wave func- 
tion method DFT code (GPAW) 25 ' 26 . 

In each case we have employed the Perdew-Zunger local 
density approximation (LDA) for the exchange correlation 
(xc)-potential 27 . An electronic temperature of kgT m 0.1 eV 
was used to converge the Kohn-Sham wavefunctions, with all 
energies extrapolated to K. The electronic density was cal- 
culated using al2xl2xl Monkhorst-Pack special fc-point 
mesh, i.e. by using 19 special points in the irreducible Bril- 
louin zone. 

In both SPPW and PWscf, LDA based pseudopotentials for 
Li, Mg and B 28 were used, and the energy spectrum was found 
to be converged with a 25 Ry plane wave cutoff. For GPAW 
we employed a grid spacing h * 0.25 A in the MgB 2 (0001) 
surface plane, and h « 0.20 A normal to the surface, which 
yielded converged results. 

As initial atomic coordinates we consider the system to be 
a hexagonal lattice with lattice parameters a and c taken from 
the bulk experimental values, a * 5.8317 a.u. m 3.086 A and 
c w 6.6216 a.u. * 3.504 A from Ref. 1, which were previously 
shown to change by less than 0.6% upon structural optimiza- 
tion of the MgB 2 bulk 10 . Structural optimization is performed 
within QE 24 , and the system is structurally optimized until a 
maximum force below 0.001 Ry/a.u. w 0.026 eV/A was ob- 
tained. 

Schematics of the relaxed MgB 2 (0001) surface supercell 
models with B, Mg, and Li termination are shown in Fig. 1. 
The B terminated surface was modeled using a supercell con- 
sisting of 9 Mg and 10 B alternating layers, while the Mg ter- 
minated surface was obtained by adding Mg layers on top of 
each surface, with the Li terminated surface formed by replac- 
ing the two Mg surface layers with two Li layers. For all ter- 
minations four unit cells of vacuum, i.e. 4c * 26.4864 a.u. w 
14.016 A, were used to separate the surfaces. This was found 
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FIG. 1: Schematic of the structurally optimized MgB2(0001) surface 
unit cells with B, Mg, and Li termination, repeated twice in each 
direction. 



TABLE I: GPAW calculated formation energies AEf for B, Mg, and 
Li terminated MgB 2 (0001) surface in eV/A 2 and J/m 2 relative to bulk 
MgB 2 and the chemical potential of the terminating metal in the bulk 
for species X, fi x , taken from Ref. 29. 



MgBj(OOOl) 


B Terminated 


Mg Terminated 


Li Terminated 


£ s iab (eV) 


-164.31 


-170.45 


-171.92 


^Buik (eV) 


-16.95 


-16.95 


-16.95 


N f 


9 


10 


10 


fix (eV/atom) 


5.81 


1.51 


1.63 


Cb. Cug, Cu 


1,0,0 


0, 1,0 


0,-1,2 


AE f (eV/A 2 ) 


-0.01 


0.03 


-0.04 


AE f (J/m 2 ) 


-0.12 


0.56 


-0.63 



to be sufficient to ensure a good description of both surface 
and bulk states. 

To estimate the relative stability of the B, Mg, and Li ter- 
minated MgB2(0001) surfaces, we consider the surface for- 
mation energy AEf, which is defined as the total energy dif- 
ference from the bulk, and the chemical potential for each of 
the terminating species, per unit area of the surface 30 . More 
precisely, 



AE f = 



£siab - NfEsulk - T*X C XMX 



(1) 



where £siab is the total energy of the slab, ~ 29.4524 a.u. 
* 8.2475 A 2 is the surface area of the unit cell, Nf is the num- 
ber of MgB2 bulk layers per unit cell, cx is the number of 



atoms added or removed from the slab to obtain a given sur- 
face termination and fix is the chemical potential for species 
X G {B,Mg,Li}, as taken from Ref. 29. Values for all param- 
eters are provided in Table I. 

To differentiate between surface states, vacuum states, sub- 
surface states, and bulk states, we first consider the wave func- 
tion's projected density in the z-direction normal to the surface 
£>n,k(z)- This is denned as 



s DP*" 



kCx.y.z) dxdy, 



(2) 



where {x,y} are surface coordinates, S is the surface area of 
the unit cell, and i//„^ is the n th Kohn- Sham wave function 
at fc-point k. We then apply the criterion that the wave func- 
tion's weight on the surface, s,,^, surpasses a certain thresh- 
old, where we define the surface region as two bulk unit cells 
below the surface layer, and one unit cell of vacuum above the 
surface layer. In more detail, we define 



rzo+c /->L-zo+2e 

Vk = I Qnx(z)dz+ I Q n .k{z)dz, 

<JZ(j-2c J L-zq-c 



(3) 



where zo is the surface layer coordinate and L is the length of 
the unit cell in the z-direction, as depicted in Fig. 1. States 
which have more than two thirds of their weight in the sur- 
face region, i.e. s,, ^ > 0.66, are then considered to be surface 
states. 

In order to further distinguish between subsurface states 
which penetrate deep into the crystal, surface resonant states, 
and bulk states, we increased the thickness of the supercell, by 
adding six bulk unit cells in the center of the slab. The result- 
ing structure was then relaxed within QE employing a coarser 
A: -point mesh of 8 x 8 x 1, to obtain the ground-state geometry. 

To understand how surface states disperse through the SBZ 
and hybridize with bulk states, we consider the overlap inte- 
grals between the Kohn-Sham wave functions at neighboring 
Appoints, separated by Ak. In real space this is defined simply 
as 



(<An,kl l An',k+Ak)'V 



-JV. 

J'V 



tM^k+AkWdr, 



(4) 



where r = {x, y, z) is the real-space coordinate, and *V is the 
volume of the unit cell. Then ifr„ k — > i/^k+Ak for n' which 
maximizes k(4" n ,kl(^n'^+Ak)'V • I n this way we may distinguish 
between band crossings and avoided crossings in the band 
structure, and trace the surface states throughout the SBZ. 



in. RESULTS 

The DFT estimates of the surface formation energy AEf 
shown in Table I are quite small (|A£y| < 0.05 eV/A 2 ), which 
suggests that each of the MgB2(0001) surface terminations 
considered should be stable experimentally. This is clear 
when one compares with the much higher surface formation 
energies obtained for other materials, such as metal oxides, 
which are known to be stable experimentally 30 . In particular, 
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FIG. 2: (Color online) Projected density of surface and subsurface states £>„,k(z) in the B terminated MgBi(OOOl) surface versus position normal 
to the surface plane z, and Kohn-Sham eigenenergies £„ k relative to the Fermi level s F . Results from plane wave PWscf calculations before 
( — ) and after ( — ) structural optimization are compared with real-space GPAW calculations for the relaxed structure (•••)• Solid and dotted 
vertical lines denote positions of the B and Mg atomic layers, respectively. Inserts show f?„_k(z) from plane wave PWscf calculations for a 
relaxed extended supercell model with six bulk unit cells added in the center of the slab. 



we find it should be thermodynamically preferred to replace 
a Mg surface termination by a Li termination, so that surface 
doping by Li is quite feasible. It should be noted that our 
choice of bulk metallic B, Mg, and Li for the chemical poten- 
tial's shown in Table I has a significant influence on the result- 
ing formation energies obtained. However, the main conclu- 
sion that each surface termination should be obtainable under 
particular experimental conditions should still hold. 

The projection of the bulk band structure in a direction per- 
pendicular to the crystal surface (F — > A) is a fingerprint 
of the particular crystal structure. The same is true for a 
MgB2(0001) surface, for which the projected band structure 
is actually quite simple. Due to the weak B-B layer over- 
lap, there are three rather narrow B cr bands (cr\, cr 2 , and <x 3 ), 
while the overlap between p z orbitals in the B layers and s or- 
bitals in Mg layers yields a wider n band. The n band is then 
separated from cr bands by two wide gaps, which are notice- 
able in Figs. 3, 6 and 9. 

Since the surface band structure of B and Mg terminated 
MgB2(0001) surfaces was previously studied in Refs. 11,13, 
we will focus here on the detailed classification of the sur- 



face and subsurface states and surface state resonances in Mg, 
B and Li terminated MgB2(0001). This is accomplished by 
combining both pseudopotential-based plane wave methods 
(SPPW, PWScf) and a real space projector augmented wave 
function method (GPAW) to test the robustness of the surface 
states, and geometry optimization to show how surface relax- 
ation can modify them. 



A. B-terminated MgB 2 (0001) surface 

To understand how movement of the atomic planes in the 
surface region influences the energy and character of the sur- 
face and subsurface states, the B terminated surface electronic 
structure has been calculated both with and without structural 
optimization. The relative changes in the (B-Mg) inter-layer 
separation compared to the bulk values, Ay, are provided in 
Table II. 

We find that the relative change in the inter-layer separa- 
tion is quite small, in agreement with previous plane wave 
calculations 10 using the PBE xc-functional 31 where only the 
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FIG. 3: (Color online) Band structure (lines) and energies of the sur- 
face and subsurface states shown in Fig. 2 (symbols) for the B ter- 
minated MgB2(0001) surface in eV relative to the Fermi level e f . 
Results from plane wave PWscf calculations before ( — ; ■) and after 
( — ; •) structural optimization are compared with real-space GPAW 
calculations for the relaxed structure (• • • ; a). 




r km r 

FIG. 4: (Color online) GPAW overlap-based band structure of re- 
laxed B terminated MgB2(0001) surface, with if/,, ^ — » >J/ n ',k+Ak for 
n' which maximizes |(t/'„,kl | A«',k+Ak)'v|- Surface states with more than 
two thirds of their weight in the first vacuum and first two bulk layers, 
SnM <; 0.66, are shown separately ( — ). 



TABLE II: Change in B terminated Mg-B inter-layer separation A,j 
between the i th and / h layers from the surface, relative to the bulk 
MgB2 experimental geometry. 



PWscf LDA° VASP PBE* 



A ,i[B-Mg] 


-4.4% 


-7.8 pm 


-2.1% 


-3.7 pm 


Ai, 2 [Mg-B] 


0.4% 


0.75 pm 


2.0% 


3.5 pm 


A 2>3 [B-Mg] 


-1.1% 


-2.0 pm 


0.9% 


1.6 pm 


A 3 , 4 [Mg-B] 


-1.7% 


-3.1 pm 


-1.8% 


-3.2 pm 


A 4 , 5 [B-Mg] 


-1.3% 


-2.6 pm 


0% 


0.0 pm 


A 5 , 6 [Mg-B] 


-1.3% 


-2.3 pm 






A 6J [B-Mg] 


-1.5% 


-2.7 pm 






A 7 , 8 [Mg-B] 


-1.3% 


-2.3 pm 






A 8-9 [B-Mg] 


-0.75% 


-1.4 pm 
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top four atomic layers were relaxed. This rigidity of the struc- 
ture is reflected in the band structure calculations shown in 
Fig. 3. We find that the band structures with and without struc- 
tural optimization do not differ significantly for the B termi- 
nated surface. 

The energies of the surface and subsurface states at the T, 
Kand M points for the structurally unoptimized and optimized 
surfaces are also shown in Fig. 3. We find that the B termi- 
nated surface has several surface and subsurface states, con- 
sisting of three localized bands which are well separated from 
their corresponding bulk bands. For these localized bands we 
plot g„ k in Fig. 2 at F, K and M points. The deepest subsur- 
face state's band (Fig. 2(a,b,c)) consists of B en orbitals, and 
is localized mainly in the first and slightly in the second B 
layer. Specifically, we find for this band s„ ^ > 0.66 through- 
out the SBZ, as shown in Fig. 4. In fact, the B cr\ surface band 
may be described semi -quantitatively by a +0.6 eV shift of the 
top of the B (j\ bulk bands. 

To differentiate between localized states, surface reso- 
nances, and bulk states, we recalculated the electronic struc- 
ture for a thicker slab with six additional MgB2 bulk layers. 
These results are shown as insets in the upper right comers of 
the Q n \Sx) plots in Fig. 2. For the larger slab we clearly see in 
Fig. 2(b) that the B cr\ state is indeed a surface state localized 
on the topmost B atomic layer. 

The second group of surface states (Fig. 2(d,h,j)) are a com- 
bination of Mg s and B p z orbitals. These sp- surface states 
are mainly localized above the surface and between B and Mg 
layers. Figure 4 shows that the sp z surface state is clearly 
recognizable as it disperses through the SBZ, although some 
hybridization occurs when crossing the narrow B cr 2 3 band. 
As with the B cr\ band, the sp z surface band is well described 
by a +0.6 eV shift of the top of the respective bulk band, in 
this case the n bands. We also find for B termination the sp z 
states decay slowly into bulk, which is also seen for the ex- 
tended slab shown as insets in Fig. 2(d,h,j). 

The third group of surface states (Fig. 2 (g,e,i)) are B <x 
states which are more than 90% localized in the topmost B 
layer. From Fig. 4 we see that the B 1x2,3 surface states may 
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TABLE III: Change in Mg terminated Mg-B inter-layer separation 
Ajj between the i* and f h layers from the surface, relative to the 
bulk MgB 2 experimental geometry. 



PWscf LDA" VASP PBE* 



Ao,i[Mg-B] 


-5.2% 


-9.2 pm 


-3.7% 


-6.5 pm 


A li2 [B-Mg] 


-0.3% 


-0.5 pm 


1.2% 


2.1 pm 


A 2 , 3 [Mg-B] 


-1.9% 


-3.3 pm 


0.2% 


0.4 pm 


A 3 , 4 [B-Mg] 


-0.7% 


-1.2 pm 


0.5% 


0.9 pm 


A 4 , 5 [Mg-B] 


-0.4% 


-0.7 pm 


-0.3% 


-0.5 pm 


A 5 , 6 [B-Mg] 


-0.9% 


-1.6 pm 






A 6 , 7 [Mg-B] 


-0.75% 


-1.3 pm 






A 7 , 8 [B-Mg] 


-0.9% 


-1.6 pm 






A 8 , 9 [Mg-B] 


-0.9% 


-1.6 pm 






A 9 ,io[B-Mg] 


-0.9% 


-1.5 pm 
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be clearly traced throughout the SBZ, at about 0.6 eV above 
their respective B cr 2 .3 bulk bands. Thus, Fig. 4 shows that for 
the B terminated MgB2(0001) surface the B o"i, sp z , and B 
surface bands are near quantitatively described by rigidly 
shifting the top of the respective bulk bands up in energy by 
0.6 eV. 



B. Mg-terminated MgB 2 (0001) surface 

Even though optimization of the crystal structure only 
slightly modified the Mg terminated surface (cf. Table III) 
and band structure (cf. Fig. 6), modifications of the surface 
electronic structure are quite radical. Density distributions 
in Fig. 5 show that for the unrelaxed surface there exist two 
types of localized states, an sp z surface state (Fig. 5(a)) and 
B cr subsurface states (Fig. 5(d,e)). On the other hand, even 
though the relaxed surface is only slightly compressed (max- 
imally -2.7% in the first layer), it induces the appearance of 
new surface states (Fig. 5(b,c)). Such states consist of B cr 
states, which are localized around B layers but penetrate quite 
deep inside the crystal. At the K point this state may even be 
considered a surface state resonance. Figure 7 clearly shows 
that these states are more properly associated with the B <T2 
and cr 3 surface states (Fig. 5 (d,e)). 

Another consequence of the surface relaxation is that the 
B cr surface states (Fig. 5(d,b,e)) are pushed upwards towards 
the surface B layers. In other words, in the relaxed crystal the 
cr surface states are completely localized in the one or two top- 
most B layers. As shown in Fig. 7, the three B cr bands have 
definite surface character throughout the SBZ, where they fol- 
low the bottom edge of the B cr bulk bands. 

We also find relaxation lowers the energy of the cr subsur- 
face states by approximately 0. 1 eV at the F and K points, and 
increases the energy of the sp z surface state from -1.86 eV to 
- 1 .78 eV at the F point, in better agreement with the experi- 
mental value of -1.6 eV. It should also be noted that the B cr 3 
surface band is more localized on the topmost B layer in the 



TABLE IV: Change in Li terminated Mg-B inter-layer separation 
Aj j between the i th and / h layers from the surface, relative to the 
bulk MgB 2 experimental geometry. 



PWscf LDA" 



A ,i[Li-B] 


-29.8% 


-52.4 pm 


A I2 [B-Mg] 


-0.6% 


-1.0 pm 


A„[Mg-B] 


-1.2% 


-2.2 pm 


A 3 , 4 [B-Mg] 


-1.5% 


-2.7 pm 


A 4 , 5 [Mg-B] 


-1.2% 


-2.1 pm 


A 5 , 6 [B-Mg] 


-1.5% 


-2.6 pm 


A 6J [Mg-B] 


-1.4% 


-2.4 pm 


A 7 , 8 [B-Mg] 


-1.9% 


-3.3 pm 


A 8 , 9 [Mg-B] 


-1.5% 


-2.6 pm 


A 9 , 10 [B-Mg] 


-1.5% 


-2.7 pm 
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real-space GPAW calculation. This is attributable to the in- 
herit difficulties plane wave based methods have in describing 
state localization, compared to real-space methods. 

In fact, from Fig. 7 we see that the bottom of all three 
B cr bands have s„k > 0.66 throughout the SBZ. This sug- 
gests that adding a layer of Mg shifts the B cr surface bands 
down in energy by about 0.9 eV relative to a the B terminated 
MgB2(0001) surface. On the other hand, we find that the sp z 
surface band is higher up in energy by about 0.9 eV relative 
to the B terminated surface. In Sec. IV we will discuss how 
this rigid shifting in energy of the surface bands may be un- 
derstood in terms of charging of the MgB2(0001) surface. 



C. Li-terminated MgB 2 (0001) surface 

To understand how surface charging, or electro-chemical 
doping might affect MgB2(0001) surface states, we have di- 
rectly replaced Mg by Li in the first atomic layer (Mg— >Li). 
In effect, this approximates the removal of one electron from 
each of the Mg terminated surfaces. In other words, replac- 
ing Mg by Li resembles a positively charged Mg terminated 
(Mg +I ) surface. This allows us to directly probe the effect 
of surface charging of the Mg terminated surface without the 
need for compensating background charges in the calculation. 
We have then performed structural optimization of the Li ter- 
minated MgB2(0001) surface to model how Li termination 
changes the surface states. From Table IV we see that the 
smaller atomic radius of Li induces a substantial relaxation, 
with the relaxed Li atomic layer moved down by 52.4 pm to- 
wards the B atomic layer. 

Both the surface charging and the crystal structure mod- 
ification strongly change the electronic structure of the sur- 
face states. Figure 9 shows that the lowest B cr x sur- 
face resonance band for the structurally unoptimized surface 
(Mg— >Li~Mg +I ), is slightly separated from the upper edge of 
the bulk B cr\ band. The density distributions of such states 
are shown in Fig. 8(a,b,c,f), where we see that B crj surface 
resonances are mainly localized in the four topmost B layers. 
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FIG. 5: (Color online) Projected density of surface and subsurface states £>„.k(z) in the Mg terminated MgB 2 (0001) surface versus position 
normal to the surface plane z, and Kohn-Sham eigenenergies e„ k relative to the Fermi level e f . Results from plane wave PWscf calculations 
before ( — ) and after ( — ) structural optimization are compared with real-space GPAW calculations for the relaxed structure (•••)• Solid and 
dotted vertical lines denote positions of the B and Mg atomic layers, respectively. Inserts show p„,k(z) from plane wave PWscf calculations for 
a structurally optimized extended supercell model with six bulk unit cells added in the center of the slab. 




FIG. 6: (Color online) Band structure (lines) and energies of the sur- 
face and subsurface states shown in Fig. 5 (symbols) for the Mg ter- 
minated MgB 2 (0001) surface in eV relative to the Fermi level e f . 
Results from plane wave PWscf calculations before ( — ; ■) and after 
( — ; •) structural optimization are compared with real-space GPAW 
calculations for the relaxed structure (■ ■ • ; a). 



FIG. 7: (Color online) GPAW overlap-based band structure of re- 
laxed Mg terminated MgBi(OOOl) surface, with i/f„ k — > i/v,k+Ak f° r 
n' which maximizes |(t/'„,kl | /V,k+Ak)'v|- Surface states with more than 
two thirds of their weight in the first vacuum and first two bulk layers, 
■Vk ~ 0.66, are shown separately ( — ). 



Structural optimization causes a slight downward shift of the 
surface resonance B cr\ bands and induces a strong interac- 
tion/hybridization with the bulk B cr\ band, so that they decay 
into the bulk states. I.e., it increases the amplitude of the prop- 
agation of the resonance state into the bulk, as can be seen at 
the T point in Fig. 8(a). For K and M points, the B <T\ sur- 



face resonances cannot even be distinguished from the normal 
bulk states and are not shown. This is seen clearly in Fig. 10, 
with the relaxed Li terminated B cr\ surface states becoming 
increasingly bulk-like between F and K, until s„ k < 0.66. 

Relaxation also causes a 0.4 eV reduction of the binding 
energy of the sp z supra-surface state band. This is expected 
because relaxation causes contraction of the topmost Li-B in- 
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FIG. 8: (Color online) Projected density of surface and subsurface states g„x(z) in the Li terminated MgB2(0001) surface versus position 
normal to the surface plane z, and Kohn-Sham eigenenergies e„ k relative to the Fermi level b f . Results from plane wave PWscf calculations 
before ( — ) and after ( — ) structural optimization are compared with real-space GPAW calculations for the relaxed structure (•••)• Solid, 
dotted, and dashed vertical lines denote positions of the B, Mg, and Li atomic layers, respectively. Inserts show £>„,k(z) from plane wave PWscf 
calculations for a structurally optimized extended supercell model with six bulk unit cells added in the center of the slab. 
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FIG. 9: (Color online) Band structure (lines) and energies of the sur- 
face and subsurface states shown in Fig. 8 (symbols) for the Li ter- 
minated MgB2(0001) surface in eV relative to the Fermi level e f . 
Results from plane wave PWscf calculations before ( — ; ■) and after 
( — ; •) structural optimization are compared with real-space GPAW 
calculations for the relaxed structure (• • • ; a). 



FIG. 10: (Color online) GPAW overlap-based band structure of re- 
laxed Li terminated MgB2(0001) surface, with i/f„_k — > (/v^+Ak for 
n' which maximizes |(t/'„,k|t/v,k+Ak)'v • Surface states with more than 
two thirds of their weight in the first vacuum and first two bulk layers, 
•Sn.k <: 0.66, are shown separately ( — ). 
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terstitial region {cf. Table IV) where the sp z surface state is 
mostly located {cf. Fig. 8). The binding energy of the sp z sur- 
face state at Y point is strongly reduced, from about 2.6 eV to 
0.1 eV, comparing with the B terminated surface. This sug- 
gests that the binding energy of the sp z surface state at the 
F point may be strongly influenced by surface contamination 
or doping of the B terminated surface, as we will discuss in 
Sec. IV. 

We also find the more localized B <x 2 and cr 3 bands are 
strongly modified by Li termination. As can be seen in Fig. 9 
for Mg— >Li~Mg +1 , the B o~2 and <x 3 surface bands are always 
slightly separated from the upper edge of the bulk B o"2 and 
cr 3 bands, with their density highly localized in the first one 
or two B layers, as seen in Fig. 8(g,e,h). These states resem- 
ble the highly localized B <x 2 and 0-3 surface states in the B 
terminated surface. However, surface relaxation shifts the lo- 
calized B o"2 and 1x3 bands down in energy. This causes a 
strong interaction between localized and bulk <x states, with 
the surface state decaying into the bulk states. This is clearly 
illustrated by the densities of these states for the relaxed Li ter- 
minated structure shown in Fig. 8(e,h). Here we see that the 
cr localized state now penetrates substantially into the crys- 
tal. By plotting the same distributions for an extended crystal 
structure (cf. insets of Fig. 8(e,h)) we see that these states are 
actually localized state resonances or even pure bulk states 
at the K point. From Fig. 10 we see that the B cr 2 and cr 3 
surface states strongly hybridize with the bulk <x states, with 
s„,k < 0.66 over much of the SBZ. 

From this we make two observations. First, removing 
charge from the Mg terminated surface would strongly affect 
the surface and subsurface states. Second, for the relaxed Li 
terminated surface all subsurface states hybridize into local- 
ized state resonances or bulk states, while the sp z surface state 
survives. We will discuss how these affects may be under- 
stood and controlled, with the goal of tuning the surface states 
of MgB2, in the following section. 

IV. DISCUSSION 

Comparing the relative positions of the B cr and sp z surface 
states for B, Mg, and Li terminations, shown in Figs. 4, 7, 
and 10, respectively, clear trends emerge. These suggest that 
changing surface termination or doping electro-chemically 
can tune both the energy and localization of the MgB2(0001) 
surface states. 

To demonstrate this point, we replot in Fig. 11 the ener- 
gies of the surface and subsurface states at the F point as a 
function of surface termination. As we move from left to 
right in Fig. 11, charge is transferred to the B atomic layer, 
as we move from neutral (B), to charging of less than -\e 
(Mg->Li~Mg +1 ), of about -le (Li), and finally, of about -2e 
(Mg). As charge is donated to the B atomic layer, both the 
degenerate B cr2 and cr 3 levels (shown in Fig. 11(a)) and B 
(j\ levels (shown in Fig. 11(c)) are filled and move down in 
energy relative to the bulk MgB2 Fermi level, in a quasi-linear 
fashion. Further, when surface states cross the bulk B cr bands, 
they hybridize, and lose much of their weight on the surface. 
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FIG. 11: Surface and subsurface state energies at the T point, e n f 
versus surface termination and doping of the MgB 2 (0001) surface 
by B, Mg— >Li~Mg +1 , Li, and Mg, in eV relative to the Fermi level 
ep, for (a) B <t 2 and tr 3 , (b) sp : , and (c) B <T\ states. Results from 
plane wave PWscf calculations before (■) and after (•) structural 
optimization are compared with real-space GPAW calculations for 
the relaxed structures (a). The B cr 2i 3 and B bulk band regions 
are shaded gray. Fits to the surface state positions ( — ) are provided 
as guides to the eye. 



This is clearly seen in Fig. 10, where for the B <t 2 and cr 3 
states s j= < 0.66, and the B cr 1 state shown in Fig. 8(a), also 
has more weight in the bulk. This strongly suggests that by 
electro-chemically doping the surface, we may tune both the 
energy and localization of the B cr surface states. 

For the higher energy B o"2 and cr 3 surface states this is 
particularly interesting, as they resemble quantum well (QW) 
states. This is because the occupying electrons effectively feel 
the potential of a QW centered around the B layer. More- 
over, electrons in such states are trapped in the topmost B 
layer but are "free" to move within the band. In this way, 
the B cr band may be considered a 2D electron gas decou- 
pled from the surrounding 3D electron gas of the bulk MgB2. 
Such a two-component electron plasma allows different types 
of plasma dispersions relations. For example, a noble acous- 
tic plasmon may be present, whose dispersion resembles that 
of the acoustic plasmons. Such a plasmon may affect (screen) 
the electron-phonon coupling. 

On the other hand, for the sp z state shown in Fig. 11(b), 
the energy of the surface state does not correlate with charg- 
ing of the B atomic layer. Instead, to understand the behavior 
of the sp z state, we must first consider its projected density 
g n f(z). For each termination we find that the sp z state is an 
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ionic/metallic bonding state between the Mg/Li and B atomic 
layers. This is clear from Figs. 2(d), 5(a), and 8(d), where 
we see that the sp z state has mostly p z character localized in 
the inter-layer region. Specifically, for B termination the sp, 
has significant weight in the bulk, while for Mg and Li ter- 
minations the sp, state has more s character, with significant 
weight in the vacuum region. 

This is also clear from the band structures shown in Figs. 4, 
7, and 10. Specifically, for B termination (cf. Fig. 4) the sp z 
state is at the top of the broad n band, having mostly p z char- 
acter, and is filled throughout most of the SBZ. For the Mg 
terminated surface (cf. Fig. 7) the sp z state has more s charac- 
ter, is separated from the bulk n bands although still following 
the 7i bands through the SBZ, is shifted up in energy, and par- 
tially emptied. When we directly replace Mg by Li to model 
a Mg +1 surface, we find that the sp z surface state has mostly 
Mg/Li s character (cf. Fig. 8(d)), is shifted up in energy be- 
tween the n and Mg s bulk states, and is almost completely 
empty. For the relaxed Li terminated surface (cf. Fig. 10) the 
sp z state is nearly emptied, and has significant s character, as 
shown in Fig. 8(d). 

Together, this suggests that by altering the termination of 
the MgB2(0001) surface, the sp, surface state is changed 
from a bulk-like state with p z character for B termination, 
to an increasingly localized and emptied surface state with 
more metallic s character from Mg termination, to Mg— >Li (or 
Mg +1 ) termination, and finally the structurally optimized Li 
termination. This is consistent with previous calculations for 
the Na terminated MgB2(0001) surface 13 , which found the sp, 
surface state between that calculated here for Mg and Mg— >Li 
terminations. In effect, by varying the surface termination 
or electro-chemically doping the surface, one can tune both 
the localization and energy range of MgB2(0001) sp, surface 
states. 

In this way, one could then modify the density or Fermi ve- 
locity of the charge carriers on the surface and consequently 
surface plasmons dispersions. This means that surface plas- 
mon characteristics may be tuned rather easily by surface 
charging or doping. 

V. CONCLUSIONS 

In this paper we have investigated MgB2(0001) surface and 
subsurface states, and the influence of surface termination (by 
B, Mg, and Li) and surface charging (replacing Mg by Li to 



model Mg +1 ) using both plane wave and real-space ab initio 
methods. Generally we find for the MgB2(0001) surface there 
exist three types of surface states which are often well sepa- 
rated from bulk bands, namely, the B a-\ surface band, the sp z 
supra-surface state band, and the B cr 2 and cr e QW bands. 

If charge is added to the B layer by either changing the sur- 
face termination or surface charging, the B cr surface bands 
shift shift down in energy. As these bands cross the bulk B cr 
bands, they lose their surface state character due to hybridiza- 
tion. In effect, by adjusting the charging of the B subsurface 
atomic layer, we may tune both the localization and energy 
range of the B cr surface states. 

On the other hand, the addition of a Mg, Mg +1 , or Li layer 
above the B surface layer gives the sp, surface state band an 
increasingly metallic s character as it is emptied. We find 
these results are quite robust, and independent of the ab initio 
methodology employed (plane wave or real-space) and relax- 
ation of the surface. 

Overall, these results suggest that by changing surface cov- 
erage and charging, one may tune both the density and Fermi 
velocity of the charge carriers at the surface. This would make 
possible the controlled tuning of the surface plasmon disper- 
sions, e.g. the slope of the acoustic surface plasmon, which 
may then influence (modify or enhance) "surface" supercon- 
ductivity. This has profound implications for both the super- 
conducting behaviour and surface or acoustic surface plasmon 
dispersion of MgB2. 

Therefore, because of the diversity of MgB2's surface elec- 
tronic structure, with combinations of various types of quasi- 
2D and 3D plasmas, it seems that the MgB2 surface should 
have very unusual dielectric and optical properties. To clarify 
these unique properties, our next step will be the investigation 
of the MgB2 surface dielectric response. 
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